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SUMMARY 


A method has been developed for contact problems which may be either frictional or 
frictionless and may involve extensive sliding between deformable bodies. It is based on an 
assumed stress hybrid approach and on an incremental variational principle for which the 
Euler's equations of the functional include the equilibrium and compatibility conditions at 
the contact surface. The tractions at an assumed contact surface are introduced as Lagrang- 
ian multipliers in the formulation. It has been concluded from the results of several example 
solutions that the extensive sliding contact between deformable bodies can be solved by the 
present method. 


INTRODUCTION 


The finite element method has been applied by many authors for solving solid mechan- 
ics problems which involve undetermined contact surfaces. They include the relatively 
simple Hertz contact problem for which there is no sliding between contact surfaces and the 
small displacement model can be used [1, 2]. They also include problems involving relative 
sliding either with friction or in frictionless conditions [3]. Existing solutions are based 
largely on the conventional assumed displacement finite element model . 

The present paper is based on an assumed stress approach and on an incremental vari- 
ational principle for which the Euler's equations of the functional include the equilibrium 
and compatibility conditions at the contact surface. An assumed contact surface is inserted 
between bodies in contact and is divided into elements. Contact tractions are independently 
assumed in terms of unknown values of such nodes of the contact elements. Thus, a finite 
element equation includes nodal displacements and nodal contact tractions as unknown. 

This paper is to present the variational principle and the corresponding finite element imple- 
mentation for this problem. 
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Coordinates before deformation 

Displacements defined within an element 

Displacements defined along element boundaries 

Tractions 

Contact tractions 

Coefficient of friction 

The whole boundary of the body 

The contact surface 

The portion of the boundary where displacements are prescribed 
The portion of the boundary where loads are applied. 

Volume 

Strains 

Stresses 

Quantity tangential to the contact surface 

Quantity normal to the contact surface 

Matrix 

Tensor 

Prescribed quantity 
Quantity pertinent to body A 
Quantity pertinent to body B 
Quantity pertinent to bodies A and B 
Quantity pertinent to element N 

Incremental quantity 
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GENERAL INCREMENTAL ASSUMED STRESS FORMULATION FOR CONTACT PROBLEM 


The requirements for contact (the conditions of contact) are as follows: 

(I) At the point of contact between two bodies tractions exerted 
on each other are the same in magnitude and are opposite in 
directions. 


(2) The normal tractions are compressive and the tangential 
tractions counteract relative movement of the bodies. 

(3) There should be no gap and no penetration of material points 
at the place of contact. 


A B 

Consider two bodies A and B shown in Fig. 1 with volumes V and V , and boundary 

a & A B A B 

surfaces, 3 V* and 3 V which are composed of portions, Sj- and S^-, and S^ and S^. 

These two bodies share a contact surface, through which they interact. The previously 

mentioned conditions of contact in incremental form are 
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and the signs in Eqs. 3 and 4 are chosen such that tangential tractions on this surface act to 
restrain the relative movement of contacting points. 

A finite element method which is based on a variational principle with relaxed con- 
tinuity requirement at interelement boundaries is defined as a hybrid model [4]. Boland and 
Pian [5] have applied an incremental assumed stress hybrid method for large deflection 

analyses of thin elastic structures. The functional IT that has been derived in the refer- 
/ me 

ence [5]! based on the Updated Lagrangian coordinate system is used as Ibase for deriving 

the functional for the present problem. 

The conditions of no gap and no overlapping on ijie place^of contact are introduced 

into the functional by means of Lagrangian multipliers, T + aT and T + aT . The 
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The Euler's equations of 3T are 
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On the contact surface S_ , in addition to Eqs. 5 and 6 
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Introducing frictional constraint on T + aT and T +aT , such that 
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and rearranging Eqs. 11, 12, 13, 14 and 15, results in Eqs. 1, 2, 3, and 4. Thus, 

it has been proved that the conditions of contact are Eulers's equations of the functional Tf C 
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These Euler equations are the strain displacement relations, the mechanical boundary 
condition, the stress equilibrium requirements along the interelement boundaries and the 
conditions of contact. Since they are only satisfied in an average sense within an increment 
they cannot be expected to satisfy these conditions in the usual sense. It is, therefore, 
necessary to consider a compatibility check, a stress equilibrium check, and a contact check. 

It is seen that Eq. 7 already has all of these built-in checks. The 5-th and 6-th 
integral terms in the functional correspond to the equilibrium check and 7-th term to the 
compatibility check. Also, the compatibility and equilibrium checks of the contact surface 
(contact check) are easily identified, the equilibrium check being the 8-th integral term 
and the compatibility check, the 9-th term. 


FINITE ELEMENT IMPLEMENTATION 


Since the aim of the present work is to solve contact problems by the finite element 
method, expressions arising from nonlinearities, not due to contact, are excluded from 
equations. A technique for solving these equations, with only contact nonlinearities, will 
be discussed here. 

Neglecting honlinearities not due to contact, the assumed stress hybrid functional, 
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The stresses a# are expressed in terms of a finite number of stress parameters, ap and 
the element boundary displacements interpolated in terms of the nodal displacements, 
and Aq. Also, coordinates are interpolated in terms of their nodal values. As a solid 
continuum is subdivided into elements, the contact surface is also discretized into finite 
number of elements referred to here as "contact elements" with "contact nodes. " The con- 
tact traction J + aT is interpolated in terms of its nodal values, t. 

^ Thus, interpolations ofjthem are: = P* a£, aT = R»a<j, U = L • q, 

AU = L,* Ac^ and T + AT = M»£. Substituting these interpolations into Eq . 16, 
and defining the following matrices 
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in which Q?, and Rp result from stress equilibrium checks on and on the contact 
q C N n 

surface S_ , and R_ , from the initial mismatch checks on S_ . 

N N C N 


Equation 17 contains three unknown vectors, the incremental stress parameters the 

incremental displacements acj, and the contact stresses t . Whereas, the vector are 
independent on the element level, aq and t are not. Thus, eliminating from Eq. 17, 
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Summing up over all elements and taking the variations of the functional with 
respect to and t, results in 
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Equation 19 represents the total assembled finite element matrix equation. 
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TECHNIQUE OF ITERATIVE SOLUTION 


Once the contacting bodies are adequately constrained such that the inverse of the 
global stiffness matrix, in Eq. 19 can be calculated, it can be used throughout the 
iteration procedure. Thus, in locating the contact surface only J5^ needs to be recomputed 
in each iteration. The global stiffness matrix K, remains constant during this process. Even 
in the case of material and/or large deflection nonlinearitles, it is possible to use a modified 
Newton-Raphson method,* hence the global stiffness matrix K may remain constant during 
this process. 

For a two-dimensional problem, before each iteration, the contact surface is a line, 
fixed in the coordinate system, but not to the contacting bodies. Such line is assumed 
known in order to perform the necessary integrals. But if, before the iteration, one has 
assumed the location of the contact surface and the positions on it that the contacting nodes 
of the bodies will occupy, then it would appear that the displacement increment can be 
specified. This is not the case. Instead, the problem is solved for the displacement incre- 
ments and if the contact surface found therefrom is not coincident with the one presupposed, 
a new contact surface location is calculated, and then an iteration can be followed. 


It has been found that the length of the contact element which yields best results is 
the same as the length of the contacting side of boundary elements of the bodies. In order 
to facilitate programming, the nodes of the contact element are chosen to be coincident 
with those of one of the contacting bodies. As a result of variations in load, the place of 
contact changes? thus a vital part of the solution is to establish a procedure for calculating 
this change. A trial and error scheme is employed because it is virtually impossible to form- 
ulate a variationalpr incipleincluding unknown variables which locate the surface of contact. 


The overall strategy for solving the contact problem is discussed here. First, an incre- 
ment in the external load or prescribed displacement is applied. Second, a contact surface 
is assumed together with the points on it through which nodes of the bodies contact each 
other. Also, the types of contact (sliding or non-sliding) at each of the above-mentioned 
points a re assumed. For the initial calculation of the first load increment, the above assump- 
tions are made simply by inspection, and for the first iteration after each new load increment, 
the converged solution of the previous load step is used. Third, all the necessary matrices 
are calculated and assembled. At the i-th iteration of the N-th load step, incremental 


displacement, aU and contact tractions ( T, + aT.) are solved from a finite element 
k k k N_i 

matrix equation* Fourth, knowing the total displacement U, at the end of the previous 

N-1 i k 

loading step N-1, the total displacements + a on the boundary, and contact 

tractions *(T^ + aT^) are checked to determine if they satisfy the conditions of contact. 


If they do not satisfy these conditions, the location of the assumed contact 
surface is modified and the procedure repeated until they do. Next, a 
convergence test is made, and if the solution is not convergent, the location of the 
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contact surface is further modified and the solution procedure repeated. 

To determine if the solution satisfies the conditions of contact, the following 
assurances are made: 

(1) That nodes of either body beyond the last contacting nodes from the 
previous iteration have not penetrated the other body. 

(2) That tractions at the contacting nodes are compressive. These normal 
tractions can be calculated by three different methods; (a) from the 
stress coefficients, (b) from the equivalent nodal forces, and (c) from 
the contact tractions, T + aT. Here the last method was used. 

(3) That the relationship between normal and tangential contact tractions. 

It + aT \ * “If + aT I 

s s n n 

is satisfied. Depending on which of the above checks, if any, is violated, one of the 
following procedures is employed to modify the assumed location of the contact surface. 


(a) If (1) is violated, the contact surface may be extended to include the 
points at which penetration has occurred. 

(b) If (2) is violated, the contact surface is reduced by excluding nodes 
at which the tractions are tensile. 

(c) If (3) is violated, sliding is allowed to occur. 


After the conditions of contact are satisfied, a test for convergence can be made by 
calculating the following quantity: 
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where aU^ is the displacement at the k-th degree of freedom. If R is less than a precribed 
quantity, say 0.01, the solution is considered as converged. 

EXAMPLE SOLUTIONS 

The finite element model and solution scheme are applied to problems of contact 
between a disk and a semi-infinite half-plane. The overall mesh pattern, the location of 
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the prescribed displacement and relevant dimensions are shown in Fig. 2-a, with the area 
immediately surrounding the contact surface shown in great detail in Fig. 2-b. The semi- 
infinite half-plane has been modeled by a finite one with overall dimensions much larger 
than those of the disk. The basic element used is four-node quadrilateral element derived 
by assuming seven 0 - parameters and linear displacement distribution along each edge. 
Five-node and six-node elements are also introduced in transition regions between coarser 
and finer meshes. Contact tractions along each contact element are approximated by linear 
interpolations. 

Non-Sliding Contact 

Problems are solved for the case with both applied loads and prescribed displacements 
at the top of the disk. In the problems, the ratio of Young's moduli are varied over a range 
from 1 to 10"^ and slightly different mesh patterns near the contact surface are used to 
accommodate’, node-to-node and node-fo-internode contacts. Loads or displacements are 
applied by fhethree increments until the length of contact surface becomes about 2.4 mm. 
For each increment, the converged solutions are reached with three or four iterations. For 
these solutions, the best results for contact tractions are obtained when calculated from 
equivalent nodal forces and are compared excellently with the Hertz solution in all cases. 

Frictionless Contact With Extensive Sliding 


The half disk and the semi-infinite half-plane are also used to demonstrate the capa- 
bility of this formul at i on to solve extensive sliding contact problems. Since the contact 
between the two bodies is frictionless, the solution is independent of the path; thus a Hertz 
solution is again available for comparison. Solutions are obtained for prescribed displace- 
ments at the top of the disk by eight increment. Stress distributions on the plane of con- 
tact, for two prescribed displacements, are plotted in Fig. 3 where zero position represents 
the point of initial contact. It is seen that the solution agrees almost exactly with that of 
the Hertz solution. It is noted that the center of symmetry of the stress distribution moves 
to the left as the half-disk slides in that direction. 

Frictional Contact With Extensive Sliding 

The same problem is again solved here with friction between the disk and the half- 
plane as an added consideration. No solution to this problem, analytic or otherwise, can 
be found; thus, results arrived at here will be justified by comparison with the results of the 
previous section and by showing that they satisfy the conditions of contact. 

The normal tractions at the contacting nodes between the disk and the half-plane for 
every 4th displacement increment are shown in Fig. 3 and are compared with those of the 
frictionless case. Because of friction, it can be seen that the displacement of the contact 
surface is retarded. That the normal tractions of the plane and disk are equal in magnitude 
and opposite in sign is also evident in the figure. This implies that the normal tractions 
satisfy a condition of contact. A condition of sliding contact requires the ratio between 
normal and tangential components of tractions to be constant and equal to the coefficient of 
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friction . They were approximately verified for all contacting nodes. Finally, the con- 
tacting surfaces of the disk and the plane are shown in Fig. 4 along with the locations of the 
nodes obtained in previous solutions. It can be seen that friction retards the movement of 
the contact surface, and as in the previous solution, through averaging over the entire con- 
tact surface, the contact condition of no separation or penetration is satisfied. 

CONCLUSIONS 

(1) An incremental variational principle and a corresponding finite element 
formulation have been made for contact problems based on an assumed stress 
hybrid method. An iterative scheme for the solution has been developed. 

(2) Successful applications of the present method for plane elasticity problems 
have been demonstrated for 

(a) Non-sliding problems with node-to-node contact 
and with node-to-internode contact, 

(b) Frictionless contact with extensive sliding, and 

(c) Frictional contact with extensive sliding. 

(3) The present method should be extended to problems Involving material 
and/or geometrical nonlinearities in addition to contact nonlinearity. 
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Figure 3.- Normal traction on contact surface. 
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